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ABSTRACT 

We present a sample of nine very high resolution cosmological simulations in the mass range of 
My- 1T — 7 x 10 11 — 4 x 10 12 M Q starting from ACDM initial conditions. Our simulations include 
primordial radiative cooling, photoionization, star formation, supernova II feedback, but exclude 
supernova driven winds and AGN feedback. We confirm our earlier results with higher resolution 
simulations and demonstrate that the simulated galaxies assemble in two phases, with the initial 
growth dominated by compact (r < r c g ) in situ star formation fueled by cold, low entropy gas streams, 
whereas the late growth is dominated by accretion of old stars formed in subunits outside the main 
galaxy (r > r c g). We find a positive correlation between the fraction of accreted stars and the final 
mass of our galaxies, with the most massive galaxies having the largest fraction of accreted stars. 
The two-phase formation mechanism naturally explains the observed downsizing, bimodality and size 
growth of the galaxy population. The mean assembly redshifts of the in situ made stars are very 
similar at Zf^ ns ~ 2.5 for all galaxies, whereas the assembly redshift of the accreted stars show a much 
larger scatter at Zf. acc ~ 0.5 — 1.5. Very high resolution simulations show that gravitational feedback 
strongly suppresses late star formation in massive galaxies contributing to the observed galaxy color 
bimodality. However, additional heating sources probably in the form of AGN and SNI feedback are 
also required to prevent late gas inflows and associated residual star formation in the more massive 
galaxies. The accretion of stellar material (dry minor mergers) is also responsible for the observed size 
growth of early- type galaxies. Consistent with their assembly histories we find that the dark matter 
fractions within the stellar half-mass radii continuously increase towards lower redshift from about 
/dm ~ 0.05 at z ~ 3 to /dm ~ 0.1 — 0.3 at z = 0. In addition, the logarithmic slope of the total 
density profile is nearly isothermal at the present-day (7' ~ 1.9 — 2.2) also in good agreement with 
recent lensing observations. Our simulations predict almost constant slopes until redshift z =1 and 
then steeper slopes of 7' ~ 3 at higher redshifts (z > 2). Finally, the input of gravitational heating 
lowers the central dark matter densities in the galaxies, but the effect is much smaller compared to 
the simulations without supernova feedback and in better agreement with other recent studies. Our 
highest resolution simulations demonstrate that the tendency of high mass galaxies to greatly reduce 
their star formation at late times follows from fundamental physical principles. 

Subject headings: cosmology: theory - galaxies: elliptical and lenticular- galaxies: formation - galax- 
ies: evolution - methods: numerical 



1. INTRODUCTION 

In the classic cold dark matter (CDM) picture of 
galaxy formation, gas falling into dark matter halos is 
shock-heated approximately to the halo virial tempera- 
ture, T vir = 10 6 (v C i rc /167kms~ 1 ) 2 K maintaining quasi- 
hydrostatic equilibrium with the dark matter componen t 
(jRees fe Ostrikerlll977t ISflkl 119771; IWhite fe Reeslll978l ). 
Then, given a density-profile for the g cooling radius 
can be calculated, inside of which the gas can radiate its 
thermal energy away, losing its pressure support, and set- 
tling into a centrifugally suppor ted disk while conservin g 
its specific angular momentum ([Fall fc Efstathioulll980l ). 
Eventually the collapsing gas reaches high enough den- 
sities to allow for star formation and a stellar disk is 
grown from the inside out. Over the years this simplified 
picture has been updated and extended into the pow- 
erful framework of se mi-analytic galaxy formation (e.g. 
White fc Frenlll991l lKauffman n~et al.lil9 93: Cole et all 
2000; De Lucia et al. 2006: Somcrville et al. 2008) that 



has been successful in reproducing many of the present- 
day properties of observed galaxies. 

A first challenge of the standard model o f galax y for- 
mation was already brought up by iBinnevi (|1977f ) , who 
used analytical models of protogalactic collapse to argue 
that under some plausible physical conditions the frac- 
tion of shock-heated gas could be small with the majority 
of the gas having temperature s below T v j r . This idea was 
revive d bv iBirnboim fc Dekell ([2003); Dc kel fc Birnboiml 
(2006), who assuming spherical symmetry showed that 
galaxy halos can only shock-heat infalling gas if the cool- 
ing rate for gas behind the shock is lower than the com- 
pression rate of the infalling gas. This criterion trans- 
lates to a critical minimum mass for halos that are able 
to shock-heat the infalling gas of M s hock ~ 10 116 Mq, 
with this mass being roughly independent of redshift and 
very simil ar to the critical mass isolat ed al r eady in the 
papers bv iRees fc Ostrikerl (fl977h and[sii! (fl977h . Ha- 
los below this mass are not massive enough to support 
stable shocks and the majority of their gas is accreted 
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cold. However, galaxies above the critical shock mass 
can also be fed with cold gas along filaments penetrating 
deep inside the hot halo. This general picture has since 
been confirmed in cosmological simulations using both 
Smoothed Particle Hydrodyna mics (SPH) (jKeres et al.l 
120051 120091 : iBrooks et I all 12009ft and Adaptive Mesh Re- 
finement fAMR) (|Ocvirk et al.ll2008ft techniques. 

Since the early days of galaxy studies it has been 
clear that the majority of galaxies come in two types, 
either as early-type elliptical galaxies or late-type spi- 
ral disk g alaxies as ma nifested in the classical Hubble 
sequence (jHubblel Il926ft . Several recent large statisti- 
cal surveys of the local galaxy population, such as the 
Sloan Digital Sky Survey (SDSS) and the two-degree 
Field (2df) sur veys, have made this bimodality even more 
robust, (e.g. IBell et al.l [200l iKauffmann eTaTl [200l 
iBaldrv et al.ll2004ft . The division between the two classes 
of galaxies is associated with a critical stellar mass of 
M cr it — 3 x 10 10 M Q . Galaxies below this critical mass 
are typically blue, star-forming disk galaxies that lie in 
the field, whereas galaxies above M cr it are dominated by 
red spheroidal systems with old stellar populations that 
predominately live in dense environments. 

Due to the inherent complexity and non-linearity of 
the galaxy formation process, numerical computer simu- 
lations have become the ideal tool for studying structure 
formation. Interestingly, most numerical studies starting 
from cosmological initial conditions have concentrated 
on the fo rmation of late-type Milky Way-like galax- 
ies ( e. g. ISommer-Larsen et al.l 120031: iRobertson et al. | 
2004 iGovernato et al.l 12004 l2007t IScannapieco et alj 



20091 lAgertz et al.l 120091: iPiontek fc Steinmetd 120111: 
Agertz et al.ll20lTft" T Much less effort has been spent on 



study ing elliptical gala xies dMeza et a l. 2003; Naab et ah! 
120071 hereafter NO 7; iFeldmann et al.l 12010ft which is 
somewhat surprising as they contain mor e than half the 
stars in the Universe (jFukugita et al.lll998ft with the most 
massive galaxies known all being spheroidal systems (e.g. 
iThomas et aT1l2005l: iKormendv et al.ll2009ft . 

However, a proposed formation mechanism of el- 
liptical galaxies that has been studied extensively 
is the binary merger scenario, in which ellipti- 
cal galaxies are formed from the merger of two 
spiral galaxies. This picture is moderately suc- 
cessful in producing the component of intermediate 
mass, fast-rotating disky elliptic als with anisotropic 
veloci ty distri bution s. (e.g. iBarnes fc Hernquistl 
19961: ISpringell l2000t iBournaud et al . 2005; Naa b et al 



2006at ICappellari et all 120071; IHopkins et al.1 12 008b; 
Hoffman et al.l l2009t iBois et al.l l201lUCappellari et al.l 
20111 : [Emscllcn iet alJ2011t iKrajnovic et alJl2011ft . How- 
ever, the binary merger scenario has problems in re- 
producing the massive old metal-rich slowly-rotating 
boxy ellipticals, which are charac terized by low elliptici- 
ties and more is o tropic velocities (iNaab fc Burkerdl2003l: 
Cox et all 120061; IHopkins et al.l l2008al; iNaab fc Ostrikerl 
20091: IBois et al.l 12010ft . The massive slowly rotating 
ellipticals are better reproduced in simulations set in 
the cosmological context, in which the galaxy is as- 
sembled from multipl e, hierarchica l mergers of star 



burst i ng subunits. ( Burkert et all 120081 INaab et al 



2009t iJohanss on et al 



2010; Fcld mann et al. 



,2009d hereafter J09; lOser et al 
20111: lOser ef^l20Ta 



erarchical, a naive reading of the model implies that the 
most massive galaxies residing in the most massive halos 
would form last. However, two recent key observational 
results have challenged this notion. Firstly, there is now 
very strong observational evidence that old, massive red 
metal-rich galaxies were already in p l ace at redshifts of 
z=2- 3 (e.g. iBrinchmann fc Ellisll2000l: I van der Wei et al.1 
12005ft . This result can also be seen as a manifestation 
of cosmic downsizing, in which the most massive galax- 
ies formed a significant proportion of their stars at high 
redshifts, compared to the lower mass systems that have 
been forming stars over the whole cosmic epoch (e.g. 
Glazebrook et al.l 12004 1 Juneau et al.ll2005t iBundv et al.l 
20061 : ICimatti et al.l |2006|) ~ Secondly, there now also 
exists observation al evidence for significant growth in 
both the size (e.g. iDaddi et al.)l2005l: iTruiillo et al.l l2006 



20071: Ivan Dokkum et al.l l2008t Ivan der Wei etal 



van Dokku m et al.ll2010ft and mass (e.g. IBell et al.l 



2009; 



2004; 



Drorv et al.l 12004 iFaber et al.l 12007ft of massive ellipti- 
cals since z=2-3 until the present-day. The observed 
z 2 ellipticals are also typically an order of magnitude 
denser than their local counterparts and would need, in 
the absence of mergers, to undergo significant secular 
and internal evolution since z ~ 2 in order to match 
the local galaxy population ( e.g. ICimatti et ail 120081 : 
Ivan Dokkum etal.ll2008l 120091 : iBezanson et al.ll2009ft . 

Theoretically, the observed bimodality in the galaxy 
distribution can be explained if star formation in halos 
above a critic al threshold mass of M ~ 10 12 M W is sup - 
pressed (e.g. iBower et all 120061: ICattaneo et alJl2006ft . 
The quenching mechanism needs to be both energetic 
enough to trigger the quenching and long-lasting enough 
to maintain the quenching over a Hubble time. In addi- 
tion to the que nching by shock-heated gas above a crit- 
ical h alo mass ([Dekel fc Birn boim 2006; Birn boim et al.1 
12007ft . potential q uenching mechanisms inc lude the feed- 
back from AGNs (|Ciotti fc Ostrikerl l2007h . gaseous ma- 
jor mer gers triggering starburst and/o r quasar activ- 
ity ( e glMihos fc Hernquistl 11994 119961: [Springel et all 
20051: iNaabeTa l. 2006a; Bo wer et alj 2006; C roton et all 
20061 : IHopkins etail 120071 : IJohansson et al.l l2009bl lah 
and gravitational quenchi n g by clumpy ac c retion 
dKhochfar fc Ostr kerl j2008t iDekel fc Birnboiml 120081 : 



iBirnboim fc Dekelll20ill J09). 

In order to explain downsizing, the stars incorporated 
in massive galaxies need to form rapidly at high redshifts. 
This process could then proceed through the merging 
of multiple starbursting subunits, with star formation 
being potentially ef ficiently fueled by cold gas accretion 
(jDekel et alJl^Q ga) . In some respect this scenario bears 
an u ncanny resemblance to the monolithic collapse pic - 
ture (|Eggen et al.lll962t lLarsor]ll975l: Iva^Albadall 1982ft . 
After the rapid initial dissipational formation process the 
star formation is terminated already at redshifts z=2- 
3. The subsequent late assembly phase of the galaxy 
then proceeds through minor dry (gas-poo r ) merging 
(e.g Ivan Dok kum 200 j INaab et al . 2006b; Bel l et al.1 
1200a iKhochfar fc S ilk 2009; Ho pkins et aHl2009l l201fjft . 
which would explain the growth in mass and size and 
the reduction in the central density (INaab et al.l 120091; 
Bezanson et al.ll2009tlvan der Wei et al.ll2009UOser et al l 



2010L I2012T 1. Furthermore, growth by dry merging does 



Since the growth of structure in the CDM model is hi- 



not violate stringent constraints on the galaxy colors and 
amount of young stars in massive ellipticals (jBell et all 
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TABLE 1 

Galaxy Properties at r < r v i r at z = 0. 



Galaxy 




M» W 




M DM W 


Zf T™ (b) 
*T ,JJ1V1 


A DM (C) 


r vir (d) 




< A' >( f ) 


fw (g) 

J not 


m, W 


rrj /n M 
' ' v lJ 1V1 




C2 


153.5 


19.8 


9.6 


124.1 


1.3 


3.57 


250 


376 


0.035 


0.98 


1.07 


8.6 


0.125 


U 


442.0 


46.7 


37.6 


357.7 


1.5 


7.92 


356 


117 


0.038 


0.99 


29 


232 


0.25 


Y 


219.0 


24.2 


17.0 


177.8 


0.8 


4.72 


281 


405 


0.056 


0.98 


46 


368 


0.25 


A2 


185.2 


22.0 


9.4 


153.8 


1.5 


3.84 


266 


383 


0.021 


0.98 


1.29 


10.3 


0.125 


Q 


90.7 


11.5 


5.9 


73.3 


1.6 


1.74 


210 


298 


0.056 


0.89 


8.5 


69 


0.25 


T 


136.7 


14.0 


7.3 


115.4 


1.5 


2.62 


241 


310 


0.050 


0.99 


8.5 


69 


0.25 


E2 


136.4 


17.1 


9.4 


109.9 


1.6 


2.94 


240 


376 


0.047 


0.97 


1.07 


8.6 


0.125 


L 


78.6 


9.1 


5.7 


63.8 


1.7 


2.16 


200 


304 


0.030 


0.93 


8.5 


69 


0.25 


M 


79.7 


11.9 


5.5 


62.3 


0.7 


2.14 


201 


305 


0.065 


0.68 


6.0 


46 


0.25 



Note. — (a) Total masses M in 10 10 Mq; (b) Redshift at which half of the final DM component was assembled; (c) Overdensity of DM 
within a sphere of radius r = 2 Mpc, pdm('' < 2Mpc)/ < pdm > (d) Virial radius in kpc, defined as the radius enclosing an overdensity of 
200 times the critical density p cr it; (e) Maximum total circular velocity in km/s; (f) Time-averaged halo spin parameter; (g) Fraction of 
hot gas (T > 2.5 X 10 5 K); (h) Particle masses m in 10 s Mq; (i) Gravitational softening lengths in physical kpc, eoM = 2xe, 



[200l iFaber etalll2007h . 

In this fir st paper (Paper I) building on previous work 
of N07; J09: lNaab et all ()2009D we explore in more detail 
the cosmological formation of massive ellipticals using 
high-resolution numerical simulations to determine if the 
two phase picture p rovides a useful interpretation (see 
also lOser et al.ll20lc[ ). We expand on our previous simu- 
lation sample by including a total of 9 galaxies simulated 
using cooling, star formation and feedback from type II 
supernovae. Confirming our earlier studies we find that 
the simulated galaxies do assemble in two phases: The 
initial growth is dominated by compact (r < r c g) in situ 
star formation fueled by cold gas flows, whereas the late 
growth is dominated by accretion of stars formed in sub- 
units outside the main galaxy. The accreted stars assem- 
ble predominantly at larger radii (r > r e s) explaining 
both the size and mass growth of the simulated galaxies 
in broad agreement with the observations. In a com- 
panion paper (Paper II Johansson et al. in prep) we 
study the photometric and kinematic properties of our 
simulated galaxies and perform a detailed comparison 
between the simulations and recent observations. 

This paper is structured as follows. In <J2] we discuss 
the numerical code, simulation physics and parameters 
employed in this study. We begin SJ3] by discussing the 
assembly history of the dark matter and stellar compo- 
nents of the simulated galaxies. We then continue by 
discussing the size evolution and dark matter fractions 
of the galaxies as a function of their assembly histories. 
In SjHwe discuss the thermal history of the gas component 
in our simulated galaxies. Here we also demonstrate the 
importance of cold flows in the early assembly history of 
the galaxies and discuss the relative importance of su- 
pernova and gravitational feedback as a heating source. 
Finally, we summarize and discuss our findings in fJSJ 

2. SIMULATIONS 

2.1. Numerical code 

The simulations presented in this paper were per- 
formed using the parallel TreeSPH-code GADGET-2 
(jSpringell [2005). The code follows the gas dynamics 
using the Lagra ngian Smoothed Particle Hydrodynam- 
ics (SPH) (e.g. lMonaghan]|1992l) technique formulated 
in such a way that energy and entro py is manifestly 
conserved (|Springel fe Hernquistl I2002T ). The code in- 
cludes star formation and radiative cooling for a pri- 
mordial mixture of hydrogen and helium, with the rates 



computed assuming that t he gas is opticall y thin and 
in ionization equilibrium (Ka tz et all Il9 96). We in- 
cluded a spatially uniform redshift-de pendent UV back- 
groun d radiation field with a modified lHaardt fe Madaul 
19961) spectrum, where reionization takes place at z ^2 6 



Dave et alii 19991 ) and the UV background field peaks at 



z ±2 2 — 3. For a detailed investigation on the effects of 
varying the background radiation field o n the evolution of 
the galaxies p resented in this paper, see lHambrick et all 
(2009. l2011af) . 

For the simulations presented in N07;J09 we turned 
off the feedback from supernovae by eliminating the two- 
phase description of the ISM of star-forming particles in 
order to better understand galactic evolution in the ab- 
sence of any energetic feedback from non-gravitational 
sources. However, in the study presented in this paper we 
in clude the full self-regu l ated s upernova feedback model 
of ISpringel fe Hernquistl (|2003l ) in all our simulations. In 
this model the ISM is treated as a two-phase medium 
(|McKee fe Ostriker|[l977l I Johansson fe Efstathioull2006ft 
in which cold clouds are embedded in a tenuous hot gas at 
pressure equilibrium. Stars form from the cold clouds in 
regions were n > nth — 0.205 cm -3 with the short-lived 
stars supplying an energy of 10 51 ergs to the surrounding 
gas by supernovae. The threshold density, n t h, is deter- 
mined self-consistently in the model by requiring that the 
equation-of-state (EOS) is continuous at the onset of star 
formation. The star formation rate in this model is set 
by dp+/dt = (1 — f3)p c /ti,, where /3 is the mass fraction of 
massive stars (> 8M ), p c is density of cold gas and t* is 
the star formation time scale set by t* = i°(n/n t h) _1 ^ 2 . 
Each gas particle can produce a maximum of two stellar 
particles per particle, with the resulting stellar particles 
thus having half the mass of the original gas particles. 
Finally, we require an over-density contrast of A > 55.7 
for the onset of star formation in order to avoid spurious 
star formation at high redshift. 



2.2. Initial conditions 

The initial conditions of the ACDM model were created 
assuming scale-invariant adiabatic fluctuations, with 
the post-recombin ation power spectrum based on the 
parametrization of E fstathiou et all (|1992l ) with a shape 
parameter of T = 0.2. Throu ghout this paper we use a 
WMAP-1 (|Spergel et al.ll2003ft cosmology with a slightly 
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TABLE 2 

Galaxy Properties at r < r ga i at z = 0. 



Galaxy 




f (b) 
Jacc v ' 


f. (b) 




A/gas (a) 


AW 


r,W 




age* ( g ) 


f (h) 
J conv 


/dm' 1 ' 




SSFR( k » Q/SFW 


C2 


14.0 


0.73 


0.27 


1.0 


0.18 


0.96 


2.05 


2.71 


11.2 


0.56 


0.11 


0.3 


-11.6 Q 


U 


25.1 


0.81 


0.19 


1.4 


1.78 


1.00 


3.11 


4.16 


11.6 


0.35 


0.20 


3.5 


-10.9 Q 


Y 


15.5 


0.70 


0.30 


1.1 


1.11 


0.87 


2.59 


3.40 


11.4 


0.44 


0.29 


2.2 


-10.8 Q 


A2 


15.1 


0.56 


0.44 


1.6 


0.17 


0.83 


2.36 


3.19 


11.3 


0.49 


0.19 


0.4 


-11.6 Q 


Q 


8.44 


0.64 


0.36 


1.9 


0.35 


0.27 


1.94 


2.53 


11.2 


0.58 


0.15 


0.6 


-11.2 Q 


T 


11.1 


0.56 


0.44 


1.6 


0.31 


0.92 


2.44 


3.23 


11.7 


0.48 


0.20 


0.3 


-11.6 Q 


E2 


11.2 


0.64 


0.34 


1.6 


0.18 


0.90 


1.49 


1.93 


11.4 


0.51 


0.08 


0.3 


-11.5 Q 


L 


7.00 


0.35 


0.65 


2.5 


0.34 


0.99 


1.57 


2.02 


11.4 


0.54 


0.18 


0.9 


-10.9 Q 


M 


8.37 


0.61 


0.39 


0.9 


0.24 


0.12 


1.69 


2.20 


10.6 


0.67 


0.12 


0.6 


-11.2 Q 



Note. — (a) Total masses M in 10 10 Mq; (b) Fraction of accreted and in situ formed stars; (c) Redshift at which half of the final stellar 
component was assembled; (d) Fraction of hot gas (T > 2.5 X 10 s K); (e) Projected stellar half-mass radius in kpc; (f) 3D stellar half-mass 
radius in kpc; (g) Average stellar age in Gyr; (h) Baryon conversion factor; (i) Fraction of dark matter within r^jy\ (j) Mean star formation 
rate in Mq/jt averaged over the last 1 Gyr; (k) Logarithm of the specific star formation rate SFR/M* in yr — 1 ; (1) Classification of galaxy 
into quiescent (Q) or star-forming (SF) based on the SSFR 
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Fig. 1. — Dark matter mass accretion histories within the galaxy virial radii (left panels), together with the evolution of the halo spin 
parameter A' (right panels) for the full galaxy sample. Major mergers can be seen as sudden increases in both the virial DM mass and the 
halo spin parameter. With the exception of halo M all galaxies have a relatively smooth accretion histories for the last 5 Gyrs. 
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lower Hubble parameter of h = 0.65 1 with 08=0.86, 
f b = n b /n DM =0.2, O =0.3, and A o =0.7. 

The galaxies presented in this paper were simulated at 
high r esolution using the y olume renormalization tech- 
nique (jKatz fc Whitdll993t N07) by selecting target ha- 
los at z — from a low-resolution (128 3 particles) dark 
matter simulation (Lbox = 50 Mpc). The target halos 
were re-simulated at high resolution by increasing the 
particle number to 100 3 and 200 3 gas and dark matter 
particles within a cubic volume at redshift z = 24 con- 
taining all particles that ended up within the virialized 
region (conservatively we assumed a fixed radius of 0.5 
Mpc) of the halos at z = 0. The tidal forces from parti- 
cles outside the high resolution cube were approximated 
by increasingly massive dark matter particles in 5 nested 
layers. The simulated high resolution region was not con- 
taminated by massive particles. A full description of the 
methods used to generate the initial conditions can be 
found in I Weil e t al.1 (|l998h . 

In N07 we performed simulations of three galaxies at 
100 3 and one galaxy at 200 3 resolution in low-density 
environments. In the present study we expand on this 
sample by including both lower- and higher-mass galaxies 
with our final sample covering a mass range of M v ; r = 
7 x 10 11 - 4 x 10 12 M Q at z = 0. Our simulation sample 
consists of a total of nine simulations, with six models run 
at 100 3 (Halos L,M,Q,T,U,Y) and three at 200 3 (Halos 
A2,C2,E2) resolution, see Table [J for details. 

2.3. Simulation parameters 

The gravitational softening length for the 100 3 runs 
was fixed in comoving units of e gas = e* = 0.25 kpc and 
£dm = 0.5 kpc until z — 9, and after this the soften- 
ing remained fixed in physical coordinates at 0.25 kpc 
for gas and stars and 0.5 kpc for dark matter. The 
gravitational softening lengths in the 200 3 simulations 
were set to half compared to the 100 3 simulations, i.e. 
e gas = e* = 0.125 kpc and 6dm = 0.25 kpc . 

Our highest (lowest) resolution runs have particle 
masses of m sta rs = 1.07 x 10 5 M© (m sta rs = 4.5 x 10 6 M Q ), 
m gas = 2.14 x 1O 5 M (m gas = 9 x 10 6 M Q ) and m dm = 
8.6 x 10 5 M Q (m dm = 3.7 x 10 7 Af Q ) for the stars, gas and 
dark matter, respectively (see Table Q}. The SPH prop- 
erties of the gas particles were averaged over the usual 
GADGET-2 spline SPH kernel using ~ 40 ± 5 SPH gas 
particles and we imposed a minimum hydrodynamical 
smoothing length equal to 0.1e gas . 

We set the parameters governing the multi-phase feed- 
back model as follows: The star formation timescale 
£ = 1.5/i- ^yr, /3 = 0.1 for a Salpeter IMF (jSalneterl 
|1955|) . the cloud evaporation parameter A = 1000 and 
the supernova 'temperature' Tsn = 10 s K that reflects 
the heating rate from a population of supernovae for a 
given IMF. These parameter choices result in a star for- 
mation rate (SFR) that is compatib le with the observed 
Kennicutt relation (Ke nnicuttlll998f ). 

The simulations were run with high force accuracy 
Qforc e = 0-005 and time integration accuracy rj acc = 0.02 
(see ISpringell 120051 for details) using the Woodhen ma- 
chine, a 768 Core Dell Beowulf cluster hosted at the 
Princeton PICSciE HPC center. The 100 3 resolution 
simulations required about ~ 5000 CPU hrs on 32 CPUs 

1 h defined such that H =100h kms _1 Mpc _1 . 



each, whereas the three 200 3 resolution simulations were 
very expensive with each requiring about ~ 175,000 
CPU hrs on 64 CPUs. 

3. FORMATION HISTORIES OF EARLY- TYPE GALAXIES 
3.1. Assembly of the dark matter component 

We begin by studying the dark matter (DM) mass ac- 
cretion histories of our galaxy sample as shown in the 
left panel of Fig. [1] where we plot the evolution of the 
virial DM masses as a function of time. The virial ra- 
dius is defined as the radius enclosing an overdensity of 
200 times the critical density p C rit- Most of the galaxies 
assemble their DM component rapidly at high redshifts 
with typical formation redshifts of the DM component 
being Z{ t ou ~ 1.5 — 1.7, where Z{ t uM is defined as the 
redshift at which half of the final DM component was 
assembled (see Tabled]). Halos U and E2 experience ma- 
jor mergers at z ~ 1.5, whereas halos A2, Q, T and L 
show strong merging activity only at very high redshifts 
of z ~ 2 — 3 followed by relatively quiescent evolution 
until the present-day, resulting in typical DM formation 
redshifts of 1.5 for these six galaxies. The remaining ha- 
los (C2, Y and M) all experience relatively late major 
mergers at (z ~ 1.2, C2; z < 0.5, Y and M) and thus 
show lower DM formation redshifts. All galaxies with 
the exception of halo M have relatively smooth DM ac- 
cretion histories for the last 5 Gyrs. 

Our simulated galaxies were selected in relatively low- 
density regions, such that the nearest massive halo with 
M > 2 X 1O 11 M0 would be at a distance larger than 
lh~ Mpc. In order to quantify the environments of our 
simulated galaxies more accurately we calculated the 
mean overdensities of DM in the neighborhood of our 
target halos. The overdensities were extracted from the 
original low-resolution simulation by comparing the to- 
tal DM density in sphere with r=2 Mpc centered on the 
halos with the mean DM density in the simulation box. 
The resulting Adm are listed in Table [1] As expected 
this quantity is correlated with the final halo mass, with 
the most massive halos (U and Y) residing in the most 
overdense regions. However, this correlation is not mono- 
tonic, halo Q which is located in the least overdense re- 
gion does not end up with the lowest final DM mass, as 
one would naively expect. 

In the right panel of Fig. [T] we plot the evolution of 
the halo spin parameter A' defined as 



V2M viv v c r vh - ' 

where v c is the circular veloc ity at the virial rad ius r v - 1T 
and M v j r is the virial mass (jBullock et alJl200l . The 
mergers experienced by halos Y and M at z < 0.5 and 
halo C2 at z ~ 1.2 are clearly identified by sudden jumps 
in their corresponding spin parameters. At high redshift 
z > 3 the evolution of the spin parameter is very spiky 
for all the halos indicating strong merging activity. The 
final spin parameter values are consistent with the me- 
dian value o f A' ~ 0.035 derive d from large numerical 
simulations ([Maccio et al.l [2007) with the exception of 
halos Y and M, which due to their relatively late merg- 
ers end up with higher spin values (in the case of halo M 
significantly higher at A^ = 0.1). In Table [T] we also give 
the time-averaged spin parameters for all the halos. The 
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Fig. 2. — The gas masses within the galaxy radii, r ga i = r v i r /10 (left panels), together with the star formation histories computed from 
stellar ages at z = (right panels) for the full galaxy sample. For all galaxies the SFR peaks at z ~ 4 — 5 after which the rate steeply 
declines. At low redshifts most galaxies have low gas masses and corresponding very low SFRs, with the exception of galaxies U and Y, 
which have SFRs of a few solar masses per year caused by late gas inflows. 



effect of mergers incr easing the spin parameter (see also 
IVitvitska et "aT1l2002D can again be seen for halos Y and 
M, but interestingly also halos with apparently similar 
accretion histories such as halos A2 and Q can end up 
having time-averaged halo spin parameters that almost 
differ by a factor of three. 

3.2. Assembly of the baryonic component 

In Fig. [3 we study the baryonic assembly of our galax- 
ies by plotting the total gas mass within the galaxy radii 
and their star formation rates as a function of time. We 
define the galaxy radius as r ga i = r v ; r /10, meaning that 
this radius increases with time as the virial radius of the 
galaxy grows. In general the gas masses of the galax- 
ies are below 1O 1O M0 during their evolution with typical 
final gas masses of < 3 x 10 9 M Q (see Table H]). No- 
table exceptions are the most massive galaxies U and 
Y, which have significant late gas inflows at z < 1 and 
end up with final gas masses in excess of 10 10 M Q . The 
majority of the gas is hot, whe re hot gas is defin ed as 
gas with T > 2.5 x 10 5 K (see iKeres etail l2005[ ) . and 
thus not directly available for star formation. Interest- 
ing exceptions are galaxies Q and M, which both contain 
significant fractions of relatively cool gas at z = 0. 



In the right panel of Fig. [5] we show the star formation 
histories of the galaxy sample computed from stellar ages 
at z — 0. All galaxies start forming their stars in bursts 
at z = 4 - 5 with peak SFRs of ~ 60 - 120 M Q /yr for 
the more massive systems and around ^30 — 50 M Q /yr 
for the least massive systems, respectively. Subsequently 
the SFRs decline almost exponentially, with this decline 
being faster for the high-mass systems compared to the 
lower-mass systems. This general trend is only inter- 
rupted by occasional mergers that show up as bumps in 
the star formation histories. At z = 1 the typical star 
formation rates range from ~ 2 — 10 M /yr and by z — 
they have declined to below ~ 1 M©/yr for all galaxies, 
but the most massive systems (U and Y) , which still ex- 
hibit residual star formation rates of a few solar masses 
per year. 

In addition, we calculate the specific star formation 
rate (SSFR) defined as the SSFR = SFR/M,, where 
SFR is the total star formation rate and M* the total 
stellar mass of the galaxy. The value at z = can be 
found in Table [2j where the current SFR has been cal- 
culat ed using the mean SFR averaged oyer th e last 1 
Gyr. iFranx et al.1 (|2008t ): IWilliams et all (|2009f ) use the 
specific star formation rate as an indicator of quiescent 
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Fig. 3. — The evolution of the total stellar mass (top panels), the in situ formed stellar mass (middle panels) and the accreted (ex situ 
formed) stellar mass within r ga i for our galaxy sample. The arrows indicate the time at which half of the corresponding stellar component 
was assembled. In all cases the in situ component is formed early, whereas the accreted component is assembled later. 



red and dead galaxies by denning these galaxies to be 
systems with SSFR < 0.3/tjr, where tu is the Hubble 
time at a given redshift. This criterion corresponds to 
log SSFR = —10.7 at z = for our chosen cosmology. 
Thus, according to this criterion all of our galaxies can 
be classified as quiescent systems, as they exhibit SSFRs 
below this threshold value at the present-day (see Table 
[2]). The very low final star formation rates are signifi- 
cant given that the simulations include neither type la 
supernova feedback nor AGN feedback. 

3.3. The two phased formation of the stellar component 

Foll owing N07 (see also lAbadi et al.ll2003t lOser et al.l 
120101 iFeldmann et al.ll2010D we study the origin of the 
stellar component in our galaxies by dividing the stel- 
lar mass into an in situ component formed within the 
galaxy (r < r ga i) and an accreted stellar component 
formed outside the galaxy (r > r ga i), where again the 
galaxy radius is defined as r ga i = r V i r /10 at the corre- 
sponding time. The galaxies in our sample are organized 
into three groups depending on whether their late accre- 
tion history (z < 2) is dominated primarily by dissipa- 
tionless minor merging (mostly accreted stars), a mixed 
disspationless/dissipational (mostly accreted and some 
in situ stars) or a primarily dissipational (significant in 
situ stars) formation history. We find that this classi- 
fication of the assembly history divides the simulated 
galaxies broadly into more massive early-type galaxies 



with strong size evolution forming dissipationlessly and 
less massive later-type galaxies with weaker size evolu- 
tion and a more dissipational formation history. The 
evolution of the accreted and in situ stellar components 
together with the total stellar mass is shown for our en- 
tire galaxy sample in Fig. [3l The galaxies are grouped 
into three groups of three galaxies depending on whether 
the majority of the final stellar mass was accreted (C2, 
U & Y), whether in situ and accreted stellar mass were 
of roughly equal importance (A2, Q & T) or whether in 
situ star formation played a dominant role (E2, L & M, 
see Table [5]) . Note that although both galaxies M and 
E2 have more accreted stellar mass at z = 0, this is in 
both cases due to a single major merger and that the 
stellar masses in both galaxies M and E2 prior to their 
respective mergers were clearly in situ dominated (see 
Fig. HJ). 

The fraction of in situ to accreted stellar mass corre- 
lates with galaxy mass, with the more massive galax- 
ies being typically dominated by accreted stars, whereas 
the least massive galaxies in our sample have significant 
in sit u formed stellar components (see also lOser et all 
2010). However, having said that significant variations 
exist for galaxies with very similar masses, relatively 
high-mass galaxies such as galaxy A2 can have signif- 
icant in situ stellar components and a low-mass galaxy 
such as galaxy Q can have a significant contribution from 
accreted stars. Thus, the final total stellar masses alone 



Johansson et al. 



Halo U 
Halo Y 
Halo A2 
Halo C2 
Halo E2 
Halo T 
Halo Q 
Halo M 
Halo L 






Y//////// 















Halo U 
Halo Y 
Halo A2 
Halo C2 
Halo E2 
Halo T 
Halo Q 
Halo M 
Halo L 





y ///////////////A 


1 , z=l 


I 






Bins 




□ acc- 



















2 3 4 5 

Mass [10'°M S ] 



1.0 1.5 2.0 

Mass [10 10 MJ 




Fig. 4. — The absolute contribution of in situ formed stellar mass (solid colors) and accreted stellar mass (dashed colors) for our galaxy 
sample shown at z = (top left), 2 = 1 (top right), 2 = 3 (bottom left) and 2 = 5 (bottom right). All galaxies assemble rapidly at high 
redshift through in situ star formation. The late evolution is dominated by accreted stars, with the more massive galaxies ending up with 
the largest fraction of accreted stars. 



do not decide the faith of the galaxies rather the details 
in the accretion history (i.e. in situ or accretion domi- 
nated) are of a significant importance in determining the 
final properties of the galaxies. 

In Fig. 2] we show the absolute contribution of the in 
situ and accreted stellar components to the total stellar 
mass at four redshifts. Using this figure and the arrows 
in Fig. [3] indicating when half of the corresponding stel- 
lar component was assembled a clear formation picture 
emerges. The stellar component in the galaxies forms in a 
two phased process. At high redshifts (z > 3) the stellar 
component in all galaxies is assembled rapidly through in 
situ star formation, with this star formation being fueled 
by cold flows (see S]4.2[) and hierarchical mergers of multi- 
ple star-bursting subunits. At later times (z < 3) in situ 
star formation is a subdominant effect with the majority 
of the stellar growth proceeding through the accretion of 
existing stellar clumps. 

The galaxies in Fig. 0] are ordered in decreasing final 
mass from top to bottom. At very high redshifts (z = 5) 
there is no clear correlation between the stellar mass at 
this redshift and the final stellar mass at z = 0. However 
by z — 3 and z = 1 at the very latest the most massive 
galaxies at these redshifts finish as the most massive sys- 
tems at z = 0. Most galaxies have typical formation red- 
shifts of their total stellar components of Z{ .* ~ 1.4 — 1.6 
similar to their DM components (see Tableland Fig. [3]), 
where zt t * is defined again as the redshift at which half of 
the final stellar mass was assembled. However, interest- 
ing exceptions exist, galaxies C2, Y and M, which have 
significant late accretion have Zf,* ~ 1 and galaxies Q and 
L, which have very rapid initial growth result in Zf * > 2. 
The formation redshifts of the in situ made stars are very 
similar at Zf^ ns ~ 2.5 for all galaxies, whereas the forma- 
tion redshift of the accreted stars show a much larger 



scatter at Zf. aC c ~ 0.5 — 1.5. Again the final stellar mass 
alone is not the decisive factor, instead there are signif- 
icant differences in the individual accretion histories for 
different galaxies in a similar mass range. 

Finally we plot the mean ages of the in situ and ac- 
creted stellar components as a function of time in Fig. [5J 
where the stellar ages are always normalized to the age 
of the Universe at the corresponding redshift. Study- 
ing this figure a clear systematic trend can be discerned. 
The mean age of the accreted stellar component is al- 
ways older than the in situ formed stars. This can be 
understood naturally in a hierarchical Universe as the 
majority of the accreted stars form individually in low- 
mass systems at high redshifts, whereas in situ star for- 
mation proceeds at a slower rate in the more massive 
systems. Thus, the accreted stars are typically formed 
earlier and are older compared to the in situ compo- 
nent. Furthermore, as discussed above the most massive 
galaxies typically have a higher fraction of accreted stars 
and hence have older mean stellar ages. The notion of a 
two phased galaxy formation process thus naturally ex- 
plains the counter-intuitive concept of downsizing. Mas- 
sive galaxies form their stellar core mass in situ and then 
accrete substantial amounts of stars that were formed 
even earlier in smaller subsystems. Hence by z ~ 2—3 the 
most massive galaxies have the oldest stellar populations 
compared to lower mass galaxies that have still signifi- 
cant in situ star formation and much lower contribution 
from accreted stars . This trend is in good agreement with 
observations (e.g. Glazebrook et al.ll2004 I Juneau et al.l 
120051: iBundv et alj|2006l ICimatti et al. 2006) and a nat- 
ural outcome of the two phased galaxy formation pro- 
cess. At the present-day the mean stellar ages are typi- 
cally high at age* > 11.2 Gyr for all our galaxies, except 
galaxy M, which due to its very late merger has a some- 
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Fig. 5. — The evolution of the mean age of the in situ formed 
(triangles) and accreted stars (circles). The stellar ages are nor- 
malized to the age of the Universe at the corresponding redshift. 
In all cases the mean age of the accreted stellar component is older 
than the age of the younger in situ formed component. 



what younger stellar population at age„ 
Table HJ. 



3.4. 



10.6 Gyr (see 



Size evolution of early-type galaxies 

We then turn our attention to the size growth of our 
simulated early-type galaxies by studying the evolution 
of their projected half-mass radii in Fig. [BJ where we plot 
separately the evolution of R\/2 for the total stellar com- 
ponent (top panel), the in situ formed component (mid- 
dle panel) and the accreted stars (bottom panel). All 
galaxies grow in size by increasing their half-mass radii 
by factors of several between 2 = 3 and the present-day. 
However, studying Fig. [B]it is clear that the in situ and 
accreted stellar components grow very differently. For 
most galaxies the size growth of the in situ component 
is very modest with i?i/2.ms growing from ~ 0.5 kpc at 
z ~ 2 to ~ 1.5 kpc at z ~ 0. The in situ stars thus form 
typically at similar small radii in the bulge components 
of all the galaxies, regardless of the galaxy mass. The 
situation for the accreted stellar component is very dif- 
ferent with strong evolution in i?i/2,acc from < 1 kpc at 



z ~ 2 to > 3 kpc at the present-day. 

In the two phased formation mechanism the inner parts 
of the galaxies are built up primarily from in situ formed 
stars, whereas the other parts contain predominantly ac- 
creted stars, which are thus responsible fo r the observed 
size growth of early-type gal axies (e.g. iTrujillo et al.1 
120061 120071 fvan Dokkum et al .120 08. 2010). The accreted 
stars are added to the galaxies primarily through dry 
minor mergers, which are very efficient in increasing the 
final size of the system and decreasing t he fin al central 
density as demonstrat ed in Na ab et all (|2009() . see also 
iBezanson et al.l (|2009| ). Thus, we expect the amount of 
size growth to correlate with the stellar accretion his- 
tories, with the systems having the largest amount of 
accreted stars also showing the strongest size growth. 
Indeed, this is seen with the sample with a dominant 
component of accreted stars (C2, U & Y) showing av- 
erage growth in size between z ~ 3 and z ~ by a 
factor of ~ 7.9, for the samples with galaxies (A2, Q 
& T) and (E2, L & M) the average growth factors are 
~ 6.6 and ~ 5.9, respectively. In the case of galaxy M 
we also see that a late major merger is not that efficient 
in increasing the size of the galaxy, with the final R\/2 
being only marginall y larger than the pre-merger half- 
mass radius (see also iNaab et al.l I2009D . Finally, as the 
typical stellar accretion history is mass dependent favor- 
ing more accreted stars over in situ stars for increasing 
galaxy mass, we expect that the size evolution will be 
mass dependent with the most massive galaxies experi- 
encing the strongest size growth. This notion has also 
been confirmed using a larger galaxy sample simula ted 
at somewhat lower resolution (|Oser etal]l20l0ll2012l) . 

Related to the size growth we plot in Fig. [7] the cir- 
cular velocity profiles, defined as t£ = GM(r)/r for our 
galaxy sample at redshifts of z = (top panel), z = 1 
(middle panel) and z = 3 (bottom panel). We also 
give in this plot the corresponding projected half-mass 
radii of the total stellar component, demonstrating ex- 
plicitly the size growth of the galaxies from z = 3 to 
the present-day. At high redshifts the galaxies are com- 
pact with their rotation curves peaking at relatively high 
values of v c ~ 450 km/s, however none of our galax- 
ies show extreme rotation curv es with peaks in e xcess 
of v r . > 750 km/s as found by iMeza et all (|2003ft and 
Uoung et all (|2009D . With increasingly lower redshifts 
the peaks of the rotation curves are systematically shifted 
to larger radii and also lowered in absolute value demon- 
strating the size growth of the parent galaxies. At the 
present-day the peak of the rotation curve is typically in 
the range of v c ~ 300—400 km/s with only the most mas- 
sive system (galaxy U) having a rotation curve peaking 
at v c ~ 450 km/s (see Table [T]). All of our systems are 
also clearly baryon dominated in their inner regions with 
the DM fraction being typically low at /dm ~ 0.1 — 0.3 
within the central 3D half-mass radii (see Table [5]) . 

Comparing our rotation curves to the ones presented 
in N07 for simulations without supernova feedback we 
see that the rotation curves in the present study peak at 
values higher by roughly ~ 100 km/s. In addition, the 
rotation curves with supernova feedback are also some- 
what more centrally peaked indicating a more massive 
central bulge component. The main reason for this differ- 
ence lies in the different star formation threshold, which 
is higher by a factor of ten in the present study. In 
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Fig. 6. — The evolution of the projected half mass radius for all stars (top panels), for the in situ formed stars (middle panel) and for 
the accreted stars (bottom panel) for our galaxy sample. The evolution of the projected half mass radius remains relatively flat for the in 
situ component. On the contrary the half mass radius of the accreted component evolves strongly resulting in a size evolution of the total 
stellar component. 



the no feedback simulations presented in N07 stars were 
very efficiently formed already in infalling subclumps. 
In the current simulations star formation is somewhat 
delayed due to the higher star formation threshold and 
the inclusion of supernova feedback. Thus the gas is re- 
quired to reach higher densities in the central regions 
of the galaxy before it can be turned into stars, thus 
explaining the slightly higher and more peaky rotation 
curves of our current simulation sample. However, out- 
side the central (r < 3 kpc) bulge-dominated region all 
rotation curves are relatively flat showing a constant ro- 
tation curve i n relatively good a g reement with obser- 
vations (e.g. iBolton et alJ I2008t ICoccato et all 120091: 
iNapolitano et al.l 120091 l201lh ! Finally, we reproduce a 
key finding of N07 who showed that increasing numeri- 
cal resolution lowers the peak of the rotation curve and 
makes it overall flatter. In the top right panel of Fig. [7] 
we plot as an comparison the rotation curve of galaxy 
E run at 100 3 resolution (El) and compare it our stan- 
dard simulation of the same initial conditions at 200 3 
resolution (E2). Increasing the resolution from 100 3 to 
200 3 lowers the peak of the rotation curve by ~ 40 km/s 
and makes the overall shape somewhat flatter in good 
agreement with the findings presented in N07. 



3.5. Evolution of the dark matter fractions 



In Fig. [8] we study the evolution of the dark matter 
properties of our simulated galaxies by showing the dark 
matter fractions inside the three-dimensional half-mass 
radius (ran) as a function of redshift. At high redshifts 
the simulated galaxies are all very strongly baryon dom- 
inated due to gas feeding in cold gas flows and the en- 
suing rapid in situ star formation. Typical dark matter 
fractions within one r^n at z ~ 4 — 5 are very low at 
/dm % 0.05. At lower redshifts of z < 3 the importance 
of in situ star formation is reduced and the majority of 
the stellar growth proceeds through stellar dry accretion. 
For most galaxies the final dark matter fractions are es- 
tablished at z ~ 1.5 — 2 after which they remain more 
or less constant. The exceptions being galaxies U and 
Y, which both have significant late gas inflows (see §3.2|) 
that also drag in DM increasing the overall dark matter 
fractions below z < 1. The final dark matter fractions 
at the present-day for our galaxy sample within one r^n 
are in the range of /dm ~ 0.1 — 0.3 (see Table[2]). 

Our simulated dark matter fractions can be compared 
to recent observations that in a self-consistent way com- 
bine the constraints from both gravitational lensing and 
stellar kinematics and derive estimates of the dark mat- 
ter fractions of early-type galaxies (e.g. iKoopmans et alJ 



[2lMlAuger et alJ200ai2010tlBarnabe et al.l2009ll20Tl 

The derived dark matter fractions from these observa- 
tional studies are sensitive to the stellar initial mass func- 
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Fig. 7. — The circular rotation curves for our galaxy sample at z = (top panel), at z = 1 (middle panel) and at z = 3 (bottom panel). In 
addition the projected half mass radii at the corresponding redshift is given. In order to demonstrate the importance of numerical resolution 
we plot in the top right panel the rotation curves for halo E at both 200 3 (E2) and 100 3 (El) resolution. All galaxies are dominated by 
luminous matter in their central parts within a few half-mass radii. 



tion (IMF) with a Chabrier IMF typically resulting in 
dark matter fractions that are larger roughly by a fac- 
tor of two compared to the dark matter fractions d erived 
using a Salpeter IMF (see e.g. lAuger et all 12010( 1 . The 
observed dark matter fractions derived using a Salpeter 
IMF seem to cluster in the range /dm ~ 0.1 — 0.3 in good 
agreement with our simu lated results dKoopmans et aF 



20061: lAuger et all 120091 and see also iGrillo k Gobat 



27Hj1 : lGrilldl201(^ Given our limited statistics we also 
see hints of a weak trend with mass, with the more mas- 
sive galaxies having on average slightly larger dark mat- 
ter fractions of /dm ~ 0.2 with the lower mass galaxies 
clustering around a value of /dm ~ 0.1. The general 
trend of larger dark matter fractions for more massive 
galaxies i s in broad agreem ent with the observational 

data (e.g. lAuger et al1l2010D 

Following iKoopmans et all (|2006l ) we fit the total den- 
sity profile of our simulated galaxies using the profile 



(2) 



PtotM =P0 — 
JO 



where po can be uniquely determined and ro can be set 
arbitrarily (we chose ro = r^o)- The only remaining free 
parameter in the density distribution is then the loga- 
rithmic density slope 7' = —d\ogptot/d\ogr, which we 
fitted between 0.1 kpc and Sr^jj. The fitting procedure is 
insensitive to a cutoff beyond several half-mass radii, but 



significantly more sensitive to the cutoff at the inner ra- 
dius. Tests showed that the chosen cutoff radii provided 
a robust and fair estimate of the central total density pro- 
files for our simulated galaxies. In Fig. [9] we show the 
evolution of the logarithmic density slope 7' as a func- 
tion of time for our simulated galaxies. At high redshifts 
all of the galaxies are baryon dominated and the cor- 
responding density profile is very steep with 7' ~ 3. At 
lower redshift the fractions of dark matter within the half 
mass radii increases and the corresponding total density 
profiles become flatter, with all galaxies approaching the 
isothermal value of 7' ~ 2 by the present-day. The mean 
density slope of all our galaxies at z = is (7') = 2.06 
with an intrinsic scatter of A ~ 0.2. This is in excellent 
agreement with the latest observational determinations, 
which found (7') = 2.078 ± 0.027 f or a sample of 73 
early- type galaxies (| Auger et al.ll2010| ). Again, given our 
rather limited sample we cannot draw strong conclusions 
concerning the mass dependence of 7'. However, we find 
that the more massive galaxies typically have larger val- 
ues of 7' ~ 2.2 at z = 0, whereas the lower mass galax- 
ies typically show 7' values below the isothermal limit 
of 7' ~ 1.9. Again, the trend of more massive galaxies 
typically having larger v alues of 7' is in ag reement with 
recent observations (e.g. lAuger et al.ll2010f ). 

3.6. Evolution of the baryon conversion factor 
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Fig. 8. — The evolution of the dark matter fraction inside the 
three dimensional half mass radius as a function of redshift for 
our simulated galaxy sample. At high redshifts the dark matter 
fractions are typically very low at /dm ~ 0.05, with the fractions 
increasing by a factor of a few to /dm ~ 0.1 — 0.3 by the present- 
day. 
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Fig. 9. — The evolution of the total logarithmic density slope 
7' = — d log ptot / d log r as a function of a redshift for our simulated 
galaxy sample. At high redshifts the galaxies are baryon dominated 
and the total density profiles are typically steep at 7' ~ 3. At lower 
redshifts the dark matter fractions in the galaxies increase and the 
density slopes flatten towards values of 7' ~ 2 at z = 0. 



We end this Section by studying the galaxy formation 
efficiency In order to do this we define the baryonic 
conversion parameter / con v as 

/conv — In 1 ' I J 

™DM,vir/ "DM Jb x ?™DM,vir 

where /& = ftb/^DM — 0-2 is the baryon fraction, m* lga i 
is the total stellar mass in the central galaxy and ?7iDM,vir 
is the total DM mass within the virial radius. In effect 
the conversion factor, / con v, thus tells us what fraction 
of the total available baryonic mass in the halo is locked 
in the central stellar component. In Fig. [TU]we show the 
evolution of the conversion factor as a function of time. 
At high redshifts (z > 2) the evolution of the galaxies 
is driven by multiple mergers resulting in a spiky evo- 
lution of the conversion factor reaching very high values 
of /conv ~ 0.8. The subsequent evolution is quieter with 
/conv remaining more or less constant below z < 1, with 
the exception being galaxy M, for which the late merger 
shows up as a marked increase in the conversion factor. 



At the present-day our galaxies have typical conver- 
sion factors of / con v ~ 0.4 — 0.6 (see Table [5]). These 
values are about a factor of ~ 2 — 3 higher than the val- 
ues predicted from halo occupation statistics, in which 
dark matter halo masses are matched to central stel- 
lar masses assuming a one-to-one mono tonic relation- 
ship (jGuo et all 12010: iMoster et alll2010D . These stud- 
ies find that the galaxy formation efficiency peaks with 
a conversion factor of / conv ~ 0.2 for halo masses of 
-Whaio ~ 6 x 10 11 M Q . Our simulated halo masses at 
z = are M halo - 8 x 10 11 - 4 x 10 12 M Q and should thus 
contain galaxies with conversion factors of / conv ~ 0.2, 
which is c learly not the case. Furth ermore weak lensing 
studies by Mandelbaum et all ()2006l ) indicate that galax- 
ies with central stellar masses of M* jCon ~ 11.2 x 10 10 A/q 
(typical stellar masses for our galaxy sample, see Table 
[2]) should have corresponding DM halos with masses of 
M D M,vir ~ 34±g° x IQ^h^MQ, which would result in 

conversion factors of / con v ~ 0.16lg;J 5 . 
The conversion factor declines rapidly both towards 
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Fig. 10. — The evolution of the baryon conversion efficiency, 
i.e. the fraction of available baryonic mass locked in stars in the 
central galaxy. Typical values of / conv ~ 0.4 — 0.6 are a factor of 
2-3 higher than expected from halo occupation statistics. 



lower and higher mass es around the p eak halo mass of 
Mhaio - 6 x 1O U M (|Guo et all 12010) . Towards lower 
masses this decline is most probably caused by UV and 
X-ray heating and supernova driven winds, which are 
able to exp el significant amounts of gas f r om dwarf galax- 
ies (e .g. lOppenheimer fc Pavel 120061 : lHambrick et al.l 
l2011al) , whereas for systems more massive than the peak 
halo mass the reduction in / conv is possibly due to an in- 
creasing contribution from feedback from supermassive 
BHs at t he centers of the galax ies, which is able to expel 
gas (e.g. iMcCarthv etalll2010l) . 

These effects are not included in the present study, 
thus it is not too surprising that our simulated conver- 
sion factors are too high by a factor of 2-3. The main 
question is how the inclusion of these effects would in- 
fluence the relative contribution of in situ formed and 
accreted stars to the final stellar component. Supernova 
driven winds are most efficient in suppressing star forma- 
tion in small systems and hence would primarily affect 
the accreted stellar component that is formed on average 
in sm aller galaxies ([Governato et al.l [20101 : lAgertz et al.l 
120111: iGuedes et al.l 1201 ll) . However, also more massive 



systems are influenced by supernova driven winds and 
thus also the formation of the in situ component would 
be slowed down. AGN feedback is only important in 
massive galaxies and would thus preferentially inhibit 
late in situ star formation in the most m a ssive galaxies 
dTevssier et al.l I20TH: iDebuhr et all I20TH; iDubois et ail 
l2012f h Thus, although the inclusion of both these effects 
would certainly lower the overall total stellar masses, the 
relative contribution of in situ and accreted stellar mass 
might not be significantly affected. The effect of dif- 
ferent feedback prescriptions on the resulting conversion 
factor is an interesting issue that deserves further study, 
but is outside the scope o f the present study (see also 
iPiontek fc Steinmet dl20Tlh . 

4. EVOLUTION OF THE GAS COMPONENT IN 
EARLY-TYPE GALAXIES 

4.1. Temperature structure of the gas 

In this Section we discuss the evolution of the gaseous 
component of early-type galaxies by using the evolution 
of halos C2, U, A2 and E2 as representative examples. 
We begin by summarizing the evolution of the gas tem- 
peratures in Fig. [TT] for the galaxies at redshifts of 
z = 0,1,3,5 by showing the temperature profiles (left 
panel), the entropy distributions (middle panel) and the 
phase-space diagrams (right panel) for all gas within the 
virial radius. The temperature of the gas is increasing 
with decreasing redshift from typical values of T ~ 10 5 K 
at z = 5 to temperatures above T ~ 10 6 K at z = 0. 
There is also a clear trend with halo mass, with more 
massive galaxies (halo U) showing higher gas tempera- 
tures indicating that the halo gas is maintained at ther- 
mal virial equilibrium. Towards the center a drop in 
temperature can be seen indicating cooling gas, however 
in the very central parts of the galaxies the temperature 
is rising again due to the input from supernova feedback 
energy in starforming particl es. 

In the middle panel of Fig. [TTJwe plot the distribution 
of the gas entropy, defined as S — kTn^ 2 / 3 for the four 
galaxies. The cooling time is defined as 



S 



10 keVcm 2 



3/2 1.5( Me //i) 2 ' (10 keVcm 2 ) 3 / 2 
(kTy/ 2 A(T,Z) 



(4) 



where k is the Boltzmann constant, A(T, Z) is the cool- 
ing function, and fj, ~ 0.59, /i e ~ 1.1 are the mean molec- 
ular weight and mean molecular weight per electron, with 
the values given for a fully ionized gas. The first term 
in the equation is a measure of the entropy and there- 
fore an adiabatic invariant, whereas the second term only 
depends on the temperature (T) an d metalicity (Z) of 
the gas through the cooling function (Sc annapieco fc Ohl 
12004 iKhalatvan et al.ll2008t ). The second factor in Eq. H 
has an absolute minimum, corresponding to t coo \ ~ 2 Gyr 
for our primordial cooling function. In Fig. [TT] we also 
plot the entropy values corresponding to minimum cool- 
ing times of 0.1, 1 and 10 Gyr as dashed lines as well as 
the entropy of neutral gas at T ~ 10 4 K marked with 
<Smax,SF- Gas with entropy below this limit is rapidly 
cooling and destined to form stars. At high redshifts the 
entropy distribution of the gas is bimodal with cold, high- 
density, star-forming gas forming a low entropy peak and 
lower density, hot shock-heated gas forming a high en- 
tropy peak. At lower redshifts the majority of the low 



14 



Johansson et al. 



E-i 
o 



Z = 



gas 





■ Halo 


C2 




■ Halo 


U 


9.6 


■ Halo 


A2 


37.6 


■ Halo 


E2 



= 9.4 
= 9.4 

100 200 300 
Radius [kpc] 




50 100 150 200 250 
Radius [kpc] 



M gas =2.1 
M gas =3.9 
M gas =3.5 
M gas =1.8 



20 40 60 
Radius [kpc] 



80 





0.4 


% 

00 




2 

\ 


0.3 










log 


0.2 


-o 




\ 

2 


0.1 







0.0 

0.4 
0.3 
0.2 

0.1 
0.0 




-4-3-2-10 1 2 3 
log S [kev cm 2 ] 




-4-3-2-10 1 2 3 
log S [kev cm 2 ] 



0.4 




-4-3-2-10 1 2 3 
log S [kev cm 2 ] 




5 10 15 20 25 30 35 
Radius [kpc] 



4-3-2-10 1 2 3 
log S [kev cm 2 ] 



7 

,— , 6 
E- 5 
4 
3 



oe 
o 









-?*W»97.5% 




f hol =99.7% 




f hol =98.3% 




f hol =97.1% 







oe 
o 



H 5 
4 
3 







,— , 6 



Ml 
O 



E- 5 
4 
3 

i 

7 
6 



ao 
o 



E- 5 
4 
3 



2 4 6 8 
log p/<p> 



10 




2 4 6 8 10 
log p/<p> 

















f hol =45.2%- 




f hol =56.8% 




f hol =61.3% 




f hol =41.2% 



2 4 6 8 10 
log p/<p> 













f hol =17.1% 
f hol =19.7% 
f hot =19.4% 
f hol =9.79% 



2 4 6 8 10 
log p/<p> 



Fig. 11. — Evolution of the temperature profile (left panel), the entropy distribution (middle panel) and the phase-space diagram (right 
panel) of all gas within r vlr as a function of rcdshift (top to bottom) for halos C2 (blue), U (orange), A2 (black), and halo E2 (red), where 
Mgas is the virial gas mass in units of 1O 1O M0 and /j lot is the fraction of diffuse gas with T > 2.5 X 10 s K and p < p t ^. The temperature 
of the diffuse gas is steadily increasing with decreasing redshift and the entropy distribution is bimodal at high redshifts. 



entropy gas has been consumed by star formation, how- 
ever in contrast to the no feedback simulations studies 
in J09 some low entropy gas capable of star formation 
remain at intermediate redshifts of z ~ 1. Having said 
that, the majority of gas (> 90%) by z ~ 1 and virtually 
all of the remaining gas by z ~ is dilute shock-heated 
gas with correspondingly long cooling times of the order 
of tcooi ~ 5 - 10 Gyr. 
In the right panel we show the phase-space diagrams 



for our galaxies, where the density is given as the density 
over the mean baryonic density at the corresponding red- 
shift. Also here the bimodality of the gas distribution can 
be seen with gas divided into a low-density hot compo- 
nent and a higher density star-forming cold component. 
The fraction of hot gas (defined as T > 2.5 x 10 5 K and 
P < Pth) is increasing steadily with decreasing redshift 
from / hot < 20% at z = 5 to / hot < 50% at z = 3. 
Starforming gas with densities above the star formation 
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threshold pth is clearly visible in the plot as most of 
the particles can be found on an equilibrium curve in 
the p — T plane dictated by the se lf-regulated feedback 
model (jSpringel fc Hernquistl l2003h . In this subresolu- 
tion model feedback from supernovae adds thermal en- 
ergy to the ISM, pressurizing the equation-of-state (EOS) 
and making it stiffer with respect to an isothermal EOS. 
The effective EOS can be written as 

PeK = (7 - l)pUoff (5) 

where 7, p and u e g are the ratio of specific heats, total gas 
density and effective specific thermal energy, respectively. 
The effective energy is the weighted average of a cold 
gas component kept fixed at T ~ 1000 K at a typical 
mass fraction of x ~ 0.9 and a hot supernova heated gas 
component at T ~ 10 s K with a typical m ass fraction 
of x ~ 0.1 (see ISpringel fc Hernquistl 120031 for further 
details). Following lRobertson et al.ll2004l we fit the EOS 
with a third order polynomial that is accurate within 1% 
for densities above the star formation density threshold 
of n t h > 0.205 cm" 3 or logn t h > -0.69, 

logPeff - 0.031 (logn H ) 3 - 0.228 (lo g 7i H ) 2 + (6) 
1.902 logriH - 11.5, 

where the pressure and number density are given in 
cgs- units, [P e e] = erg cm -3 and [n#] = cm" 3 , respec- 
tiv ely. Differences i n the n umerical factors with respect 
to iRobertson et al.l (|2004f ) are due to the fact that we 
use slightly different values for the cooling rates of a pri- 
mordial gas adopted from iTheuns et al.l (|1998| ) , which 
result in a somewhat different threshold density and cor- 
responding fit to the stiffened EOS. A few gas particles 



can be seen lying above this equilibrium curve. These 
are particles that have just received a significant energy 
injection from supernova feedback and are now cooling 
onto the equilibrium curve on the temperature decay 
scale set by the self-regulated feedba ck parameters (see 
Eq. 12 of ISpringel fc Hernquistll2003ft . 

Finally, we can derive estimates of the amount of 
cold molecular gas in our galaxies indirectly by us- 
ing the crude approximation that about ~ 90% of the 
star-forming gas above the density threshold of n t h > 
0.205 cm" 3 is in a very cold phase (T ~ 1000 K) with 
the remaining ~ 10% being supernova heated gas as 
dictated by our equilibrium supernova feedback model 
(jSpringel k, Hernquistl 12003). At high redshifts the gas 
dynamics is driven by cold gas flows and the correspond- 
ing star formation fractions are very high at /sf ~ 
50 — 60% (see ij4.3[) . As a result the predicted masses of 
cold molecular gas at z — 5 are of the order of a few times 
10 9 M Q . At the present-day the starforming gas fractions 
for all of our simulated galaxies are below < 1% and as 
result the expected cold molecular masses at z — are 
of the order a few times 10 Mg in good agreement with 
recent observations of th e typical molecular masses of 
local early- type galaxies (|Young et al.ll20lit ISerra et all 

unni). 

4.2. The role of cold and hot gas flows 

We then turn our attention to studying how the gas 
is accre ted onto the galaxies. Following the approach 
taken in lKeres et"all (|2005ll2009l ) we study the maximum 
temperature attained by accreted gas particles in order 
to deduce if they were accreted in a cold or hot mode. 
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Fig. 13. — Evolution of gas surface density, S gaa (left panel), the mass-weighted mean gas temperature, T (middle panel) and the 
mass-weighted mean gas entropy, S = kTn~ 2 ^ (right panel) for the high resolution Halo A2 at redshifts 2 = 0, z = 1,2 = 3, 2 = 5 from 
top to bottom. The length scale of each panel is depicted by the white bar in kpc, with each panel representing a slice with depth 100 kpc 
centered at the position of the central stellar component and binned with 256 2 pixels. 



In Fig. [12] we show the maximum temperatures of all 
gas particles accreted since the previous snapshot, which 
corresponds to a physical timestep of At ~ 300 Myr at 
z = 3, z = 5 and At ~ 700 Myr at z = 1, z = 0. We 
record the maximum temperature of the accreted gas 
particles at all previous snapshots excluding the snap- 
shots when the gas particles were starforming, in which 
case the high gas temperature was due to supernova feed- 
back and not hydrodynamical processes. We then relate 



these maximum temperatures to the virial temperature 
of the halos at the time of accretion and plot the relative 
mass accretion rates as a function of T max /T v i r in Fig. 

Ha 

Inspecting the figure one can see that virtually all of 
the accreted gas particles at z = 5 have temperatures 
below a tenth of the virial temperatures of the halos. 
Keeping in mind that the virial temperature at z ~ 5 
is typically a few times 10 5 K this indicates that the 
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rate dominates over the cooling rate at all redshifts for the low-density non-starforming gas. 



gas is accreted in a cold phase with typical temperatures 
of a few times 10 4 K. At z = 3 the temperature dis- 
tribution is more bimodal with the majority of the gas 
being accreted cold, but with a significant component 
of hot gas being accreted with T ~ T vil . This coin- 
cides with an increase of hot halo gas and halo masses 
of M ha i = 3 - 5 x 10 n M Q (see Fig. [T}. By z ~ 1 
the majority of the gas is accreted in the hot phase with 
some residual cold accretion and by z ~ the accretion 
of cold gas has all but disappeared. Thus the transi- 
tion from a cold gas accretion mode to a hot gas accre- 
tion mode occurs at z <~ 2 — 3, corresponding to halo 
masses of M halo = 5x lQ^M p, — l0 12 Mp, in good agree- 
ment wi th the predictions of [Dekel fc BirnboinJ (|2006f ). 
see also iRees fc Ostrikerl (|1977j ) and lSilkl (|1977| ). Below 
these redshifts the galaxy halos are massive enough to 
support stable shocks and most of the accreted gas is 
shock-heated close to the virial temperature of the ha- 
los, whereas at higher redshifts the halos are not massive 
enough to support stable accretion shocks and most of 
the gas is accreted cold. This transition is also mass- 
dependent, with an earlier transition occurring for more 
massive halos. Thus at z ~ 3 the more massive halos 
U, A2 have slightly larger fractions of hot gas than the 
lower mass halos C2, E2 (Fig ITT]) , with also some differ- 
ences in the hot gas fractions persisting until z ~ 1. At 
z = the hot gas fractions are very high for all halos at 
/hot > 97%, as most of the original gas has either formed 
stars or been shock-heated to very high temperatures. 

To better illustrate the importance of cold/hot gas ac- 
cretion we show in Fig. [13] the gas surface densities, 



together with the mass-weighted temperatures and en- 
tropies for the very high resolution halo A2 as a func- 
tion of redshift. The images are generated by projecting 
the properties of the SPH particles through a physical 
slice of 100 kpc alo ng the line-qf-sight using the gather 
approximation (see iDolag et al.l 120051) and binning the 
data over 256 2 pixels. The corresponding pixel resolu- 
tions are 0.39 kpc, 1.17 kpc, 2.44 kpc and 3.91 kpc for 
redshifts z = 5,3,1 and z = 0, respectively. The image 
resolution were chosen in order to resolve the virial radii 
drawn as solid circles at each redshift. At z — 5 the cen- 
tral galaxy is sitting at the intersection of gas filaments 
through which cold high-density gas is being fed into the 
galaxy (particularly well visible in the entropy plot as 
blue filaments). By z = 3 and especially z = 1 the 
amount of hot gas has increased significantly, however 
at both redshifts cold gas filaments are able to penetrate 
into the halo and feed the central g aseous structure wh ere 
star formation takes place (see also lDekel et al.f2009allbT ) . 
The accretion of cold gas is clumpy producing turbulent 
eddies in the gas and at z = 1 even an elongated ex- 
tended gaseous structure. By contrast the distribution 
of gas at z — is much smoother with almost the entire 
gas component found in a hot diffuse state with only very 
minor traces of cold clumpy gas remaining. 

4.3. Gravitational heating of the gas component 

Next we study the effect of gravitational feedback, 
which comes in many forms. Accreted lumps and streams 
of cold gas eventually come to rest due to drag forces de- 
positing their potential energy frictionally (i.e. through 
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decaying turbulence, see e.g. iBirnboim fc Dekell2011l) . 
Supersonic collisions of infalling gas with the ambient gas 
lead to propagating shock waves that deposit ener gy and 
entropy throughout the system (|Rvu et al.1 12003). And 
finally infalling satellite systems captured through dy- 
namical friction cause gaseous wakes from which ene rgy 
is transferred to the surrounding gas (|0strikerlll999l ). 

Following the analysis presented in J09 we plot the 
cooling and heating rates for diffuse g as, i.e. non- 
starforming gas with (p < pthresh) in Fig. Q31 The rates 
are estimated using the entropy e quation in GADGET- 
2 wh ich is for a given particle i (|Springel k, Hernquistl 
[2002]) . 
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where the first term depicts the external radiative cool- 
ing (or heating) of the gas and the second term gives the 
generation of entropy by artificial viscosity in shocks and 
where the internal energy is defined as u = A/ (7— . 
The infalling clumps have typically velocities of v ~ 
100 — 400 km/s, with the larger velocities corresponding 
to the galaxies at lower redshifts with larger potential 
wells. The corresponding sound speeds are a few fac- 
tors lower, typically c s ~ 20 — 150 km/s, thus giving rise 
to weak shocks with Mach numbers of 2-5. The gravita- 
tional feedback energy is then released in the weak shocks 
into the gas through dissipation of turbulence. 

From Fig. [14] we see that a substantial fraction of the 
gas (~ 2/3 at z = 5) and (~ 1/3 at z = 3) is starforming, 
and at z ~ 1 the starforming gas fraction remains at ~ 
10%. Comparing this to the no feedback case presented 
in J09 we see that the starforming fractions at z = 3 — 5 
are similar, whereas the feedback case shows starforming 
fractions higher by a factor of ~ 3 at z = 1 compared to 
the no feedback simulations. By z — the starforming 
fraction is very low at < 1%. 

The shock-induced heating rates of the diffuse gas is 
larger than the cooling rates at all redshifts, thus explain- 
ing the systematic increase in the diffuse gas temperature 
with decreasing redshift (Fig. [TTj) . The absolute values of 
the heating rates in the present feedback simulations for 
a given halo are somewhat lower than the corresponding 
rates in the no feedback simulation, whereas the cooling 
rates are very similar in the two simulation samples. This 
indicates that the effect of gravitational heating in the 
simulations with supernova feedback is somewhat lower, 
as expected, due to the overall suppression of low-mass 
systems that could give rise to gravitational heating when 
accreted onto the main galaxy. According to analytical 
estimates presented in J09 the efficiency of gravitational 
feedback scales as (AE) glav oc v%, where v c is the circular 
velocity of the galaxy. Thus gravitational heating should 
be stronger in more massive systems, which can also be 
seen in the z = panel of Fig. [TJJ the more massive halos 
U and A2 show heating rates above those of the some- 
what lower mass halos C2 and E2. An additional factor 
which needs to be taken into account when comparing 
halo U to halos A2, C2 and E2 is the fact that the lat- 
ter simulations have higher numerical resolution and are 
expected to show larger gravitational heating rates. At 
higher resolution smaller mass clumps can be better re- 
solved and thus the gravitational heating rates at higher 
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Fig. 15. — Cumulative DM mass distribution within r < 2 kpc 
as a function of lookback time (redshift) for halos C2,U,A2 and E2. 
Typically the central DM mass peaks at z ~ 2 — 3 after which it 
somewhat declines until z = 0. 

resolution are systematically higher as was shown in the 
resolution study presented in J09. 

An interesting case is the most massive halo U, which 
at z = show very high heating rates in the outer parts 
of the halo, with the heating rate dropping by an order of 
magnitude towards the central parts of the galaxy. This 
galaxy possesses a very extended hot gaseous halo and 
our results imply that most of the gravitational energy 
is deposited at large radii r > 100 kpc outside the cen- 
tral galaxy. Thus, for the most massive galaxies (U,Y) 
in our sample gravitational heating alone is not strong 
enough to overcome the cooling in the central parts of the 
galaxies at late times (z < 0.5) at which time most of the 
accretion process is over. Consequently these galaxies ex- 
perience strong late gas inflows as seen in Fig. [2j which 
also results in some residual star formation at z = (see 
Table [2j). The natural candidate for inhibiting these late 
gas inflows and residual star formati on is the feedback 
from a supermassiv e black hole (e.g. ICiotti et al. 2009; 
lOstriker et aT1l2010D . Thus, we conclude that although 
gravitational feedback is important in keeping diffuse gas 
hot, especially in the outer parts of the galactic haloes, 
some form of additional feedback, most probably AGN 
feedback, is required to stop late central star formation 
in the more massive galaxies. 

4.4. Gravitational heating of the DM component 

The typical dark matter fractions within the virial radii 
in our simulated galaxies are ~ 5/6 indicating at the 
same time that the majority of the available baryons are 
retained within the virial radius of our simulated galaxies 
and that dark matter is the dominant mass component 
within the virial radius. The gravitational energy input 
scales with the mass fraction of each component, thus in 
addition to heating the baryonic gas component we ex- 
pect that a fraction ft m / (f2 m +Ofe) = 5/6 of the available 
feedback energy goes into heating the DM component of 
the galaxies causing it to expand outwards and reducing 
its density within the galaxy. Calculations which assume 
adiabatic contraction of the dark matter component and 
ignore the effect of gravitational heating tend to overes- 
timate the central dark matter densities. 
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In Fig. [I5]we quantify this effect by plotting the central 
DM mass within r < 2 kpc in our four galaxies (circles) 
as a function of lookback time. Initially we see an in- 
crease in the central DM due to adiabatic contraction 
(jBlumenthal et al.lll98l IGnedin et aT1l2004l ) for all the 
galaxies with a peak central DM mass reached at z ~ 3. 
After the peak DM mass has been reached the central 
DM mass steadily declines due to the gravitational heat- 
ing of infalling clumps. The relative change in the central 
DM mass between the peak at z ~ 3 and the final value 
at z = is proportional to the amount of gravitational 
heating and thus to the amount of accreted stars. Thus 
this effect should be strongest for halos U and C2, fol- 
lowed by A2 and E2 in decreasing order. Indeed, we find 
that the central DM masses of the halos between the 
peak at z ~ 3 and z = are reduced by 38% for halo 
C2, 35% for halo A2 and 33% for halo E2 The picture is 
more complicated for halo U, the DM in this system first 
experiences adiabatic contraction at z ~ 3, followed by a 
gravitational heating induced expansion phase where the 
central DM is reduced and then at late times a second 
adiabatic contraction phase induced by the late inflow 
of gas and associated star formation at z < 0.5, thus ex- 
plaining the late upturn in the central DM mass. We note 
that the overall relative expansions of the central dark 
matter between z = 3 and z — are lower by factors 
of ~ 2 when compared to the no feedback simulations 
in J09, for which the central DM was typically reduced 
by 60% between z ~ 3 and z = 0. This again demon- 
strates the reduced importance of gravitational feedback 
in simulations that include supernova feedback and most 
importantly that gravitational feedback is now not able 
to completely negate the effect of adiabatic contraction. 
The final central DM mass of halo A2 is higher by a fac- 
tor of two compared to the corresponding final central 
DM mass of the DM only simulation presented in J09. 
Thus, in all our simulations including supernova feed- 
back we see the effects of adiabatic contraction in better 
agreement with the recent study of IGnedin et al.l (|2011h 

The ability of gravitational feedback to lower the cen- 
tral DM masses has important implications for the re- 
cent estimates of c entral dark matter c ontent of ellip- 
tical galaxies (e.g. I Thomas et al.l 1200 9') . Traditionally 
there has been some tension in the standard CDM model 
between numerical s imulations that pr edict cuspy dark 
matter profiles (e.g. iMoore et al.lll999l and the observa- 
tions that indicat e that many galaxi es have finite dark 
matter cores (e .g. iGentile et aLl2004D. As demonstrated 
here (see also lEl-Zant et al.l 120*01 r iRomano-Dfaz et al.l 
l2008l[2009| ) the gravitational energy release from infalling 
clumps might help in transforming an initially cuspy 
DM profile into a more cored profile, although additional 
sources not included in this study , such as strong wind- 
drive n supernova feedback (e .g. iPontzen fc Governatol 
120111 : iGovernato et~all l2012f ) and AGN feedback are 
probably also required to prevent late gas inflows and 
the associated excessive adiabatic contraction. 

5. DISCUSSION 

In this first paper of our two paper series we have ex- 
plored in detail the assembly histories of massive ellipti- 
cals by performing a suite of high resolution simulations 
starting from cosmological initial conditions. Our sim- 
ulations include primordial cooling, star formation and 



feedback from type II supernovae, but exclude supernova 
driven winds and AGN feedback. However, we stress 
that our simplified treatment does include, automati- 
cally, the important mechanism of gravitational heating 
which is often neglected in semi-analytic treatments but 
is a quite important aspect of maintaining overall energy 
conservation. Thus, the purpose of this paper was to see 
what aspects of the problems of the formation of mas- 
sive galaxies could be addressed without invoking major 
feedback mechanisms and to see if a careful and high 
resolution treatment from cosmological initial conditions 
would capture the essential aspects of the formation pro- 
cess. In a companion paper (Paper II Johansson et al. in 
prep) we will study the photometric and kinematic prop- 
erties of our simulated galaxies and compare the results 
with recent observations of early-type galaxies. 

In agreement with our previous studies (N07; J09; 
IQser et alJl201Qll20ll we find that there are two phases 
in the formation process of our simulated galaxies. The 
initial growth is dominated by compact (r < r e ff) in 
situ star formation fueled by cold gas flow resulting 
in early star formation peaking at very high values of 
- 50 - 120 M /yr at redshifts of z ~ 4 - 5. The later 
formation history below redshifts of z < 2 is dominated 
by dissipationless accretion of stars formed in subunits 
outside the main galaxy and the corresponding star for- 
mation rates are much lower with typical SFR values of 
< 1 M /yr. The accreted stars assemble predominantly 
at larger radii (r > r e g) explaining both the size and 
mass growth of the simulated galaxies in broad agree- 
ment with the observations. 

We also find a relation between the mass of the system 
and the fraction of accreted stars in the final galaxy, the 
larger the mass of the final galaxy the higher the fraction 
of accreted stars. This being the general picture, differ- 
ences between the simulated galaxies exist even within a 
similar mass range caused by variations in the individual 
accretion histories of the various galaxies. We also find 
that gravitational heating due to infalling baryonic lumps 
of stars and gas is a major physical process that typically 
releases of the order of ~ 10 60 ergs (J09) of which about 
1/6 (Clb/(Cl m + fib) — 1/6) goes to the gas and 5/6 to 
the dark matter, heating the former and pushing out the 
latter. The combined effect of gas exhaustion due to 
early efficient star formation and the energy release from 
gravitational feedback terminates star formation in most 
galaxies by z — 2, resulting in dead and red ellipticals by 
z = 1 in good agreement with the observations. 

The three main recent observational results concerning 
early-type galaxies are firstly the observed bimodality in 
the local galaxy population that postulates that galax- 
ies below a critical stellar mass of M cr it ~ 3 x lO lo A/ 
are typically blue, star-forming disk galaxies that lie 
in the field, whereas galaxies above M cr it are domi- 
nated by red spheroidal systems with old stellar pop- 
ulatio ns that predominately l ive in dense enviro nments 
(e.g. iKauffmann et all 120031 : fBaldrv et alJl2004 l. Sec- 
ondly, there is now very strong observational evidence 
that old, massive red metal-rich galaxies were already 
at place at redshifts of z=2-3 and that the most mas- 
sive galaxies formed a signi ficant proportion of their 
stars at high redshifts (e.g. iBrinchmann fc Ellisl 120001 : 
Ivan der Wei et aT] 120051 . Thirdly, recent observations 
have revealed growth in both the size (e.g. iTrujillo et al.1 
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2007 



2004; 



van Dokkum et af1l2008t ) and mass (e.g. iBell et ah! 
Faber et alJ^OOTf Tof massive ellipticals since z=2-3 



until the present-day. 

The two-phase formation for elliptical galaxies pre- 
sented in this paper naturally explains all of these ob- 
servations. In our simulations the central parts of the 
galaxies form very rapidly at high redshifts through cold 
accretion driven in situ star formation. The more mas- 
sive galaxies which are located in higher density regions 
assemble first, the later formation history of these mas- 
sive galaxies is then dominated by accretion of stars that 
were formed in smaller subunits at even earlier times. 
Thus, the more massive a galaxy is the earlier it formed 
and the older its average stellar population is due to 
the larger fraction of accreted stars thus explaining the 
downsizing observations. In our scenario the observed 
bimodality can be explained by the combined effect of 
gas exhaustion due to early star formation and gravita- 
tional feedback. In more massive galaxies the fraction of 
accreted stars is higher leading to stronger gravitational 
feedback heating and resulting in stronger suppression of 
star formation. Less massive galaxies have more ongoing 
in situ star formation and thus remain bluer throughout. 

Gravitational feedback is thus a contributing factor 
in explaining the galaxy bimodality, but also additional 
heating sources probably in the form of strong super- 
nova driven winds and AGN feedback are required to 
prevent late gas inflows and associated residual star for- 
mation. Additionally, the observed size growth of ellip- 
ticals is quite naturally explained in our model, in which 
dry minor mergers are responsible for the size growth, 
with the most massive galaxies showing the strongest 
growth in size as th ey have the largest fraction of ac- 
creted stars (see also lNaab et~aT1 120091: lOser et al.ll2010l 
2012). Finally, we are able match reasonably well the ob- 
servational constraints on the central dark matter frac- 
tions of early-type galaxies at the present-day. In the 
two-phased formation picture the forming galaxies are 
strongly baryon dominated at high redshifts with corre- 
spondingly low dark matter fractions and very steep total 
density profiles (7' = — d log p tot /d log r). At lower red- 
shifts the central dark matter fractions increase, although 
the total dark matter masses within r < 2 kpc are at the 
same time somewhat reduced due to the input of gravita- 
tional heating. At the present-day our simulated galaxies 
have central dark matter fractions of /dm ~ 0.1 — 0.3, 
with the logarithmic slope of the total density profile be- 
ing in the range 7' ~ 1.9 — 2.2 in good ag r eement with 
recent obse rvations (jKoopmans et al.ll2006l : lAuger et al.l 
l2tMl2010h . 

We also want to stress the importance of numerical 
resolution in galaxy formation studies. When compar- 
ing our simulation runs at 200 3 resolution (A2,C2,E2) 
to their corresponding runs at 100 3 resolution (A,C,E), 
we see some systematic effects. The final stellar masses 
in the higher resolution runs are typically lower by ~ 
10 — 20% and as a result the rotation curves for the 200 3 
runs peak at values that are ^ 20 — 40 km/s lower than 
their corresponding 100 3 runs with the overall shape also 
being flatter (see Fig [7]). More importantly the final star 
formation rates of all the 200 3 resolution are very low 
(SFR ~ 0.3 - 0.4 M /yr), with the corresponding 100 3 
resolution simulations having star formation rates that 
are typically larger by ~ 50%. Resolution has a more 



minor effect on the baryonic conversion efficiency, as the 
somewhat lower stellar masses in the high-resolution sim- 
ulations also correspondingly reside in somewhat lower 
mass DM halos. The derived differences are mainly due 
to the increased input of gravitational heating with in- 
creasing resolution, but comparing to the no feedback 
simulations discussed in N07;J09 we see that the sen- 
sitivity to numerical resolution, although still there, is 
significantly reduced in simulations that include super- 
nova feedback. Thus in conclusion, it is important to 
perform the cosmological simulations at both sufficient 
mas s (< 10 5 M(d) and spatial resolution (~ 0.1 k pc) (see 
also iGovernato et al.l l2010t iGuedes et all 1201 ll for very 
high resolution simulations of lower mass disk galaxies). 

Despite the many successes of our simple model sev- 
eral problems still remain. The most pressing issue is 
the fact that the fractions of available baryons in the 
halo that is converted into stars in the central galaxy 
are too high by a factor of ~ 2 compared to the ex- 
pectation from ha l o occ upation statistics dGuo et al.1 
l2010t iMoster et al.l I2010D and to som e recent simula- 
tions with stronger super nova feedback (jGovernato et al.l 
l2010t IGuedes et al.ll201lT ) . In addition the final gas frac- 
tions and resulting star formation rates while generally 
low are still somewhat too large, especially for the more 
massive galaxies in our simulation sample. Physical ef- 
fects that would need to be added include metal-line 
cooling, which would tend to make the problems men- 
tioned above worse. But more importantly we are miss- 
ing several important heating sources including AGN 
feedback from broad absorption line wi nd regions and 
the luminous AGN radiative out put Ce.g. ISazonoy et all 
2004L Uohansson et al.l 2009b; Booth fc Schavd I2009T 



Schaye et all 120101 : iDebuhr et al.1 l2010t lOstriker et al I 
2010D . radiative driven winds from the coupling 
with the luminous outp u t of young stars (e.g . 



Qppenheimer fc Pavel l2006t iPetkova fc Springell 120091 : 



Hopkins et alj|201ll) and feedback from type la super 



novae (e.g. Scannapieco et al.l 120051 : ICiotti fc Ostrikerl 



l2007tlOppenheimer fc Davell2008l T in addition to the sim- 
plified prescription of type II supernova feedback cur- 
rently included in the simulations. 

The effect of supernova driven winds and AGN feed- 
back is differential with respect to the masses of the 
galaxies, the former is most important for low-mass 
galaxies, whereas the latter being only important in more 
massive galaxies. We expect that the problem with the 
too high baryonic conversion efficiency will be mainly 
addressed by adding stronger supernova driven winds, 
whereas the problem of residual star formation in mas- 
sive galaxies is best addressed by adding AGN feed- 
back and type I supernova that are both mostly effective 
in the cores of massive galaxies. An additional factor 
for suppressing late star formation would be simply in- 
creasing the numerical resolution of the simulation as 
clearly demonstrated in NO 7. We have completed pre- 
liminary work incl uding strong supernov a wind feedback 
dGenel et al.ll2012| ) and AGN feedback (|Hambrick et al.l 
l201ibft and find that the addition of both of these ef- 
fects can plausibly solve the problems encountered in the 
present simulations. Thus, although the inclusion of all 
these feedback effects would certainly lower the overall 
total stellar masses in accordance with the observations, 
the relative contribution of in situ and accreted stellar 



Forming Early- Type Galaxies in ACDM simulations 



21 



mass might not be significantly affected and thus the 
general picture of a two-phased formation mechanism 
presented in this paper would remain valid. The next 
test of our cosmological two-phased formation scenario 
would be the inclusion of the aforementioned additional 
feedback sources in very high resolution zoom-in simu- 
lations of galaxies spanning a much wider range in both 
mass and environmental density, with the ultimate goal 
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